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ON DYNAMIC ALGORITHMS FOR FACTORIZATION 
INVARIANTS IN NUMERICAL MONOIDS 


THOMAS BARRON, CHRISTOPHER O’NEILL, AND ROBERTO PELAYO 

Abstract. Studying the factorization theory of numerical monoids relies on un¬ 
derstanding several important factorization invariants, including length sets, delta 
sets, and w-primality. While progress in this field has been accelerated by the use of 
computer algebra systems, many existing algorithms are computationally infeasible 
for numerical monoids with several irreducible elements. In this paper, we present 
dynamic algorithms for the factorization set, length set, delta set, and w-primality 
in numerical monoids and demonstrate that these algorithms give significant im¬ 
provements in runtime and memory usage. In describing our dynamic approach to 
computing w-primality, we extend the usual definition of this invariant to the quo¬ 
tient group of the monoid and show that several useful results naturally extend to 
this broader setting. 


1. Introduction 

Numerical monoids (co-finite, additive submonoids of N) and their non-unique factor¬ 
ization theory have emerged as an active area of study in recent years mmmmm- 
The primary objects in these investigations are factorization invariants, which measure 
in different ways the plural factorizations of an element n into irreducibles in the nu¬ 
merical monoid S. The set Z(n) of factorizations of a single element n e S can uniquely 
determine the monoid structure of S. As such, Z (n) is often too cumbersome to com¬ 
pute for large values of n in numerical monoids with several irreducible elements. Most 
of the major factorization invariants of interest, including the length set L(n), delta set 
A(n), ca-primality cu(n), and catenary degree c(n), are typically only computed after 
first computing Z (n). 

Discovering many of the results regarding these factorization invariants has re¬ 
lied on effectively utilizing computer algebra packages, such as the GAP package 
numericalsgps [8j. Of particular interest is finding closed forms for and analyzing 
the asymptotic behavior of these invariants, which frequently requires computing the 
invariant for numerous elements of S. As each of these individual computations gener¬ 
ally requires first computing Z(n) and then passing to the invariant of interest, finding 
patterns for numerical monoids with large numbers of irreducible elements is often 
computationally infeasible. 

Recent investigations in numerical monoids show that their factorization invariants 
tend to have eventual quasi-polynomial behavior. For example, cu-primality na and 
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maximal and minimal elements of length sets [2j have quasilinear behavior, while delta 
sets [7] and catenary degree [5]J have periodic (i.e., quasi-constant) behavior. The 
inductive nature of these results motivates the major results for this paper: several 
factorization invariants, including Z(n), L(n), A(n), and co(n), can be computed dy¬ 
namically. Moreover, the dynamic algorithms we present for L(n),A(n), and oj(ri) 
(Algorithms 3.7, 3.9, and 5.6, respectively) do not require computing Z(n), which 


significantly improves runtimes and memory usage. Additionally, as closed forms and 
asymptotic behavior is of interest, a dynamic algorithm, which computes several values 
quickly in succession, is a fruitful approach. 

In Section [2] we provide all necessary definitions, including those for numerical 
monoids (Definition |2.3[ ) and their length set and delta set invariants (Definition |2.6[ ), 
and state a recent periodicity result of interest. In Section [3j we give a dynamic algo¬ 
rithm to compute the delta set of a numerical monoid (Algorithm |3.9[ ), and compare 
this algorithm to existing algorithms [TO]. In Section [5j we give a dynamic algorithm 
to compute cu-primality in numerical monoids (Algorithm 5.6| ), after developing the 
necessary theory in Section [4] In Section [6j we prove Theorem 6.2 which relates 
factorizations and bullets (Definition 4.6), and derive an improved bound on start of 
quasilinear behavior of the cu-primality function (Theorem 6.12). Finally, Section [7] 
states several open questions. 

At the time of writing, Algorithms |3.9[ [373] and 543 are currently implemented in the 
GAP package numericalsgps [8]. All included benchmarks use this software, including 
comparisons to existing algorithms. 


Acknowledgements. The authors would like to thank Pedro Garcia-Sanchez for nu¬ 
merous insightful conversations. 


2. Background 

We begin by defining numerical monoids (Definition |2.3[ ) and introducing three of the 
main factorization invariants of interest: factorization sets, length sets, and delta sets 
(Definitions 2.5 and 2.6). In this section, and in the majority of this paper, monoids 
are written additively, and N denotes the set of non-negative integers. 


Definition 2.1. A commutative monoid M (written additively) is cancellative if for 
any a,b,c G M , a + b = a + c implies that b — c. An element u G M is said to be 
irreducible (or an atom ) if whenever u = a + b for a, b G M, then either a or b is a unit 
in M. The monoid M is said to be atomic if for every m G M, there exist irreducible 
elements u±, U 2 , ■ ■ ■, u r G M such that m = u\ + + • • • + u r . 


Remark 2.2. Unless otherwise stated, all monoids in this paper are assumed to be 
cancellative, commutative, and atomic. 

This paper will focus on the factorization theory of additive submonoids of N. 
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Definition 2.3. For relatively prime positive integers ni, 712 ,..., n k G N, the numerical 
monoid S generated by {ni, n 2 ,... , n k } is the additive submonoid of N given by 

S = (ni, n 2 ,... , n k ) = {cpni + c 2 n 2 H-h c k n k : q 6 N} C N. 

The Frobenius number F(S) = max(N \ S ) is the largest integer lying outside of S', 
and an element n G Z \ S is a pseudo-Frobenius number if n + rq G S for all n t . 

Remark 2.4. Every co-finite, additive submonoid of N is a numerical monoid gen¬ 
erated by a set of relatively prime generators {ni, n 2 ,... , n k }. In fact, there exists a 
unique collection of generators that is minimal with respect to set-theoretic inclusion. 
For this generating set, the irreducible elements of the numerical monoid coincide pre¬ 
cisely with the generators. Thus, every numerical monoid is cancellative, reduced, and 
finitely generated. Unless otherwise stated, when we write S = (ni, n 2 ,... , n k ), we 
assume that the n t constitute a minimal generating set and that n\ < n 2 < ■ ■ ■ < n k . 

Of particular interest in the study of numerical monoids is its non-unique factoriza¬ 
tion theory. We establish here notation for factorizations of elements in these monoids. 

Definition 2.5. Let S = (ni,n 2 ,... , n k ) be a numerical monoid. A factorization of 
x £ S is a /c-tuple a = (oq, a 2 ,..., a k ) € N k such that x = ayni + a 2 n 2 + • • • + a k n k , 
and the set of factorizations of x in S is given by 

Z s(x) = {a G : a \ n\ + a 2 n 2 + • • • + a k n k = x}. 

When the monoid is understood, this is often denoted Z(x). The length of a factoriza¬ 
tion a is given by |a| = Y^i=i a »- 

In a numerical monoid S, the set of factorizations Z(x) is usually cumbersome, 
especially for larger values of x £ S. Therefore, many factorization invariants are 
computed using only the lengths of factorizations. This leads us to the concepts of the 
length set and delta set. 

Definition 2.6. Given a numerical monoid S = (rq, n 2 ,... , n k ) and x G S, the length 
set of x is given by 

L(x) = {|a| : a G Z(x)}. 

If we order the elements of L(a;) in increasing order L(x) = {l\ < l 2 < ■ ■ ■ < l r }, the 
delta set of x is given by 

A(x) = {h — k -1 : 2 < i < r}. 

The delta set of S is given as the union of delta sets of all non-identity elements: 

A(S) = IJ A(x). 

xeS\{o} 
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FIGURE 1. A plot showing the delta sets of elements in the numerical 
monoid S = (6, 9, 20) from Example 4.5 Here, a dot is placed at the 
point (n,d) whenever d G A(n). 


Example 2.7. Consider the numerical monoid M = (6,9,20), which has Frobenius 
number F(M ) = 43. The element 60 G M has factorization set 

Z(60) = {(0, 0, 3), (1, 6,0), (4,4, 0), (7, 2, 0), (10, 0, 0)}. 

This produces the length set L(60) = {3, 7, 8, 9,10} and delta set A(60) = (1,4}. 


While the delta set of a numerical monoid S is the union of delta sets of all of its 
(infinitely many) non-identity elements, one needs only to compute A (a:) for finitely 
many x G S. A weaker version of Theorem 2.8, which appeared in [7\, states that delta 


sets in a numerical monoid are eventually periodic and that A (S') is obtained by taking 
the union of A (n) over every n < 2 kri 2 n\ + n\n *.. More recently, the authors of [ID] 
give a greatly improved bound Ng (significantly smaller than the bound from [7]) for 
the start of this periodic behavior, and show that A (S') can be computed by taking 
the union of A(n) over every n < Ns + rik — 1. 


Theorem 2.8. [TUI Corollary 18] Fix a numerical monoid S' = (ni,... ,nk). There 
exists an integer Ns such that n > Ns implies A(n) = A(n + lcm(n 1 , n^)). 

Example 2.9. Figure [I] plots the delta sets of elements in M = (6,9,20). The bound 
given in [TU] for the start of periodic behavior is Ns = 144, and we can see from the 
plot that the actual start of this behavior is 91. Additionally, the eventual period is 20, 
which by ms is guaranteed to divide (but clearly need not equal to) lcm(6, 20) = 60. 
Using this, we can compute the delta set of M to be A(M) = (1, 2, 3,4}. 


3. Factorization sets and delta sets in numerical monoids 

This section provides dynamic algorithms to compute factorization sets, length sets, 
and delta sets for numerical monoids (Algorithms |3.3[ [377] , and 3.9 respectively), along 
with corresponding proofs of correctness. Furthermore, we demonstrate the significant 
runtime advantages of these dynamic algorithms when compared to existing (non¬ 
dynamic) algorithms. 
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We begin with a lemma describing an inductive relationship between the factorization 
sets of elements. As this result applies to a significantly larger class of monoids than 
numerical monoids, it is stated in full generality. In what follows, e, denotes the ith 
unit vector in 


Lemma 3.1. Fix a reduced, finitely generated monoid M (written additively) with 
irreducible elements mi,m 2 , ■ ■ ■ , m*,. For each non-zero x G M, we have 

Z(x) = Utu {a + e ?; : a G Z(x - m,)} 

= Utr { a + e i '■ a ^ Z(x — mi), aj = 0 for each j < i}. 


Proof. Fix a G Z(x), and let i = min{j \ a j 7 ^ 0}. Then a — e, G Z {x — mf). Moreover, 
if a — ej G Z(x — mf) for some j, then j > i by minimality of i. This proves the second 
equality, from which the first equality follows. □ 


Example 3.2. Let M = (6,9,20). Lemma 3.1 allows us to compute Z(60) in terms of 


Z(40), Z(51), and Z(54). In particular, each factorization of 54 yields a factorization 
of 60 with one additional copy of the irreducible 6 . Similarly, each factorization of 51 
or 40 yields a factorization for 60 with one additional copy of 9 or 20, respectively. We 
give the full computation below. 

Z(40) = {( 0 , 0 , 2 )}, 

Z(51) = {(1,5,0),(4,3,0),(7,1,0)}, 

Z(54) = {(0,6,0), (3,4,0), ( 6 , 2,0), (9, 0,0)}, 

Z(60) = {(0,0,3)} U {(1,6,0), (4,4,0), (7, 2,0)} 

U {(1,6,0), (4,4,0), (7, 2,0), (10, 0,0)} 

= {(0, 0,3), (1, 6 , 0), (4,4, 0), (7, 2,0), (10, 0,0)}. 


Since each factorization of 60 has at least one irreducible, the above computation 
produces the entire set Z(60). Given n > 0, Algorithm |3.3| uses this method inductively 
to compute Z{i) for all i G {0,1,..., n}, with Z(0) = {0} as its starting point. 


Algorithm 3.3. Given n G S = (rq,. .., nf), computes Z(m) for all m G [0, n] D S. 
function FactorizationsUpToElement(S', n) 

F[0] 4- {0} 

for all m G [0, n] D S do 

Z ^ {} 

for all i = 1, 2,..., k with m — Hi G S do 
ZG7U{a + e,:aG F[m — rq]} 

end for 
F[m\ Z 
end for 
return F 
end function 
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s 

n 

|Z(n)| 

Z (n) 

{(to, Z(m)) : m < n} 

Alg. 

3.3 


(10,17,19,25,31) 

1000 

20293 

641 ms 

55342 ms 

20120 ms 

(51,53,55,117) 

5000 

1299 

54 ms 

61643 ms 

3874 ms 


(7,15,17,18,20) 

1000 

75375 

1028 ms 

234339 ms 

102857 ms 

(100,121,142,163,284) 

30000 

16569 

2437 ms 

13266033 ms 

660320 ms 


Table 1. Runtimes for computing Z (n) and {(m, Z (m)) : m <n} given 
G S. The last column gives runtimes for Algorithm 3.3| implemented 


n 


in GAP, and the two previous columns use the implementation of Z (n) 
found in the GAP package numericalsgps [8]. 


Remark 3.4. Table [T| gives a runtime comparison for Algorithm |3.3| with the exist¬ 
ing implementation of Z (n) in the GAP package numericalsgps. Notice that Algo¬ 
rithm 3.3| is slower than the current GAP implementition for computing a single factor¬ 
ization set Z(n), but faster for computing a large collection of factorization sets, such as 
{(m, Z(m)) : m < n}. The same holds true for memory usage: Algorithm |3.3| consumes 
more memory than the GAP implementation when computing a single factorization set, 
since it must store several factorization sets, but when computing {(m, Z(m)) : m < n}, 
it does not consume any more memory than is necessary for the returned list. 


Typically, a length set for an element is computed only once its factorization set is 
known. Lemma |3.5 provides a dynamic algorithm that allows the computation of the 
length set of an element simply by knowing the lengths sets of smaller elements (and 
not their factorization sets). 


Lemma 3.5. Fix a reduced, finitely generated monoid M with irreducible elements 
mi, m 2 ,..., mk- We have 

k k 

L(x) = |J (L(t - rrii) + 1) = {l + 1 : l G L(x - m,)} 

i =1 Z—1 

for each nonzero x G M. 


Proof. Apply the length function to both sides of the first equality of Lemma IO D 


Example 3.6. We resume notation from Example 3.2 Lemma 3.5 specializes the 


computation of Z (n) in Lemma 3.1 to length sets, allowing us to compute L(60) from 
L(40), L(51), and L(54). The key observation is that producing a factorization for 60 
from a factorization of 40, 51 or 54 using Lemma 3T always increases length by exactly 
one. Given below is the full computation. 


L(60) = (L(40) + 1) U (L(51) + 1) U (L(54) + 1) 
= {3} U {7,8,9} U {7,8,9,10} 

= {3,7,8,9,10} 
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Given n > 0, Algorithm 3.7 uses this method to dynamically compute L(m) for all 
G [0, n] fl S, which Algorithm 3.9 then uses to compute (J” =1 A(m). Notice that 


m 


this does not require computing any factorizations; see Remark 3.10 


Algorithm 3.7. Given n G S = (ni, ... ,nk), computes L(m) for all m G [0,n] fl S. 
function LengthSetsUpToElement(S', n) 

£[ 0 ] 4 - { 0 } 

for all m G [0, n] fl S do 

L 

for all i — 1, 2,..., k with rn — n t G S do 
L t— L U {7 + 1 : Z G C[m — n, ; ]} 

end for 
£[m] L 
end for 
return C 
end function 


Remark 3.8. Since Algorithm 3.7 computes {L(m 


m G [0, n] fl S}, it can also be 
m G [0, n] fl S}. As stated 
is periodic for m greater than an integer Ns described in [THlj . 


used to compute 1J^ =1 A(m) by first computing (A(m) 
A (m 


in Theorem 2.8 

Together with Algorithm |3.7 this yields Algorithm 3.9 for computing A (S). 


Algorithm 3.9. Given a numerical monoid S, computes A (S'), 
function DeltaSet(S) 

Compute C using Algorithm 8.7| 

Compute Ns using m Section 3] 

Compute A = Ume( 0 ) IV s +lcm(n 1 ,n fc )]nS A ( m ) 

return A 
end function 


Remark 3.10. In addition to giving an improved bound Ns on the start of periodic 
behavior of A s for any numerical monoid S = {n±,... , nk), the authors also give [lOl 
Algorithm 21] to find A(S) by first computing Z (Ns + rik — ni ),.... Z (Ns + rik — 1). 
Table [2] gives a runtime comparison between Algorithm 3.9 and [TOl Algorithm 21], 
demonstrating a marked improvement in computation time. 

Our method of computing A(S) has several advantages over [TUJ Algorithm 21]. 
First and foremost is memory comsumption. For large n, factorization sets grow large 
very quickly PI. whereas |L(n)| grows linearly in n p, Theorem 4.3]. Since [TOl Al¬ 
gorithm 21] requires computing Z(n) for several large n, it is often memory intensive. 
Algorithm 3.9, on the other hand, avoids the computation of factorization sets alto¬ 


gether by only computing length sets. 

Due in part to the low memory footprint, our algorithm is significantly faster than 
[ ID1 Algorithm 21]. Some runtimes for [TUI Algorithm 21] are omitted from Table [2j 


as 
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s 

Ns HD] 

Dis. 

A (S) 

HB Alg. 21] 

Alg. 

3.9 


(10,17,19,25,31) 

1180 

76 

{1,2,3} 

3254 ms 

45 ms 

(51,53,55,117) 

9699 

1006 

{2,4,6} 

23565 ms 

250 ms 

(7,15,17,18,20) 

1935 

46 

{1,2,3} 

88831 ms 

146 ms 

(7,19,20,25,29) 

3894 

76 

{1,2, 3, 5} 

— ms 

624 ms 

(11,53,73,87) 

14381 

873 

{2,4,6,8,10,22} 

49418 ms 

2588 ms 

(31,73,77,87,91) 

31364 

558 

{2,4,6} 

— ms 

4274 ms 

(100,121,142,163,284) 

24850 

5499 

{21} 

— ms 

3697 ms 

(1001,1211,1421,1631,2841) 

2063141 

114535 

{10,20,30} 

— ms 

116371 ms 


Table 2. Runtime comparison for computing A(5). Some computa¬ 
tions for m Algorithm 21] could not be completed due to insufficient 
available memory and hence have been omitted. All computations were 
completed using GAP and the package nuraericalsgps [ 8 ]. 


the high memory requirements left us unable to complete the computation. However, 
pm Table 2] gives a runtime for the computation of A((31, 73, 77, 87, 91)) on the order 
of 24,000 seconds, a stark contrast to the 4.2 seconds required for Algorithm |3.9[ 


Remark 3.11. Although Algorithm |3.9 as stated requires the computation of L(i) for 
all i < N s , Algorithm 3.7| com putes each L(i) using values at least i — n^. As such, in 
implementing Algorithm |3.9| to compute A (S'), one only needs to store rik length sets 
at any given time. The implcmenation of Algorithm |3.9 in GAP stores the length sets 
in a ring buffer of length rik, cutting the memory requirements even further. 


Remark 3.12. Since the implementation of Algorithm |3.9| in the numericalsgps pack¬ 
age i, other promising delta set algorthms have been developed. Recent results of 
the second author of this manuscript [T4j have produced an algorithm that computes 
the delta set of any affine monoid (a strictly more general setting) and appears to run 
faster than Algorithm |3.9| This algorithm is currently being implemented in GAP and 
will likely be included in a future version of the numericalsgps package. Addition¬ 
ally, a log-time algorithm to compute A(S) is given in [ 12 ] for the special case where 
S = ( 711 , 712 , 113 ) and S is non-symmetric. 


4. CU-PRIMALITY IN THE QUOTIENT GROUP q (M) 

In the remaining sections, we turn our attention from factorization sets and delta 
sets to (u-primality (Definition 4.2), a factorization invariant that has received much 
attention in recent investigations panansum. The results presented in these sections 
are motivated by those in Section [3j namely that dynamic computations can yield 
significant performance improvements when running a large collection of computations. 

In this section, we provide the necessary definitions and motivation related to the 
(u-primality invariant. We begin with a definition of the quotient group of a monoid. 
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Remark 4.1. Let q (M) denote the quotient group of a monoid M. The map M —)■ 
q(M) given by m t—)■ (m — 0) is an injective monoid homomorphism, and we often 
identify elements of M with their image in q(M) and view M C q (M). Under this 
convention, x \ y for elements x, y G q (M) if y = a + x for some a G M. 

The extension of the (n-function on a monoid to elements of its quotient group is one 
of the key insights of this paper. In this section, we develop the theory behind this 
extension, with minimal assumptions on the underlying monoid. 

Definition 4.2. Fix a monoid M. The u-primality function u>m ■ q(M) -iNU {oo} 
is given by wm(x) = m if m is the smallest positive integer with the property that 
whenever Y^i =i a * — x *= M for r > m and a* G M, there exists a subset T C {1,..., r} 
with |Tj < m such that Yhi& T a i ~ x *= Af. If no such m exists, define u>m(x) = oo. 
When M is clear from context, we simply write u(x). 


Remark 4.3. It is easy to check that the function ujm defined above coincides with 
that dehned in [ TEI Dehnitions 2.3 and 3.4] for elements of M, so Definition |4.2| simply 
extends the domain of % from M to q(M). Example 4.5 demonstrates that this is a 


natural extension, and the results that follow show that many of the properties of the 
usual (n-function still hold in this new setting. 


Remark 4.4. Written in a multiplicative setting, Definition |4.2| becomes more trans¬ 
parent; see [16l Definition 2.3] for more detail. In particular, is clear that prime monoids 
elements x G M are precisely those satisfying uj(x) = 1. In fact, much of the semi¬ 
nal work using cj-primality focused on computing u -values for irreducible elements, as 
non-unique factorizations arise from irreducible elements that are not prime. 


Example 4.5. Let S = (6,9,20) denote the numerical monoid from Example 3.2 


Since S has finite complement in N, q(S') is naturally isomorphic to Z, with the obvious 
inclusion map. Figure [2] plots side-by-side the u -values of S (as dehned in [16]) and 
those of q(S') (from Definition 4.2). Notice the plotted values coincide for each n G S, 


and the cos -values in the right-hand plot of elements lying in the complement of S seem 
to “fill in” the missing values in the left-hand plot. 


One of the key ideas used to study w-primality in recent years is its characterization 
in terms of bullet lengths m Proposition 2.10]. We now extend the definition of bullets 
to elements of the quotient group (Definition 4.6[) and recover the usual characterization 


of cn-primality in terms of their length (Proposition 4.9). Note that Proposition 4.9 
also implies that we may assume the elements a* in Dehntion 4.2 are irreducible. 


Definition 4.6. Fix a monoid M. A bullet for x G q (M) is an expression u\ + ■ • • + u r 
of irreducible elements U\,... ,u r G M such that (i) U\ + • • • + u r — x G M, and (ii) 
Mi + • • • + u r — x — Ui ^ M for each i < r. The value of a bullet u\ + • • • + u r is the 
element u\A -b u r G M, and its length is r. The set of bullets of x is denoted bul(x). 
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FIGURE 2. A plot of the u -values of elements in the numerical monoid 
S = (6, 9, 20) (left) and its quotient group (right) discussed in Exam¬ 
ple |4.5[ The smaller red dots in the right-hand plot mark ca-values of 
elements lying outside of S. 


Remark 4.7. For a monoid element m G M, the set bul(m) in Definition 4.6 is 
identical to the usual definition ([16J Definition 2.8]), since any element of q(M) that 
m divides also must lie in M. Both of these definitions differ slightly from the classical 
definition of a bullet; see 0. 

Remark 4.8. Fix a monoid M and an element x G M. If, additionally, M is both 
reduced (that is, M has no non-identity units) and finitely generated, then M has 
only finitely many irreducible elements Ui,...,Uk- In this setting, we can denote a 
bullet biU\ + • • • + bkUk G bul(x) by b = (&i, ... ,bk) G N k . This is of particular use in 
Sections [5] and |6] (as well as Example 4.11 [ ) , where ca-primality for numerical monoids 
is examined in more detail. 


Proposition 4.9. Given any monoid M , 

c o(x) = sup{r : U\ + • • • + u r G bul(a;),Uj irreducible} 
for each element x G q(M). 


Proof. This is identical to the proof of [16[ Proposition 2.10]. □ 

The notion of cover maps between bullet sets, first introduced in [15 . Definition 3.4], 
will be critical in the dynamic algorithm for computing ca-primality (Algorithm |5.6[ ). 
In the remainder of this section, we extend the original definition of cover maps to 
elements of the quotient group (Definition 4.13) and provide stronger results than the 
prior setting would allow (see Remark 4.14). First, we give Lemma 4.10, which was 
stated in [1| for numerical monoids and played a crucial role in developing the first 
algorithm to compute cu-primality. 


Lemma 4.10. Fix a monoid M, an element x G q (M), a bullet u\ + • • • + u r G bul(rc), 
and an irreducible element u G M. Then u + u\ + • • • + u r ^ bul(x). 
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Proof. Omitting u yields the expression u\ + • • • + u r , which x divides. □ 


Much like Lemma [3d" for factorizations, the bullet set of an element x is determined 
by the bullet sets of its divisors. Definition 4.13 specifies how to determine the image of 


each bullet as demonstrated in Example 4.11, and Theorem 4.12 ensures the resulting 
map is well-defined. First, we give an example. 


Example 4.11. Let S = (6,9,20) denote the numerical monoid from Example 4.5 
Consider n = 60 and the following bullet sets: 

bul(40) = {(0,0, 2), (4,4,0), (7, 2,0), (10, 0,0), (1,6,0), (0,8,0)}, 

bul(51) = {(0,7,0), (10, 0,0), (4, 3,0), (1,5,0), (0,0, 3), (7,1,0)}, 

bul(54) = {(9, 0,0), (6, 2,0), (0,6,0), (3,4,0), (0,0, 3)}, 

bul(60) = {(4,4,0), (7, 2,0), (10, 0,0), (1,6,0), (0,8,0), (0,0, 3)} 

We see that for each bullet b G bul(54), either b G bul(60) or b + ex G bul(60). Notice 
that it is impossible for both of these to lie in bul(60) by Lemma 4.10 Similarly, 
for each bullet b G bul(51), either b G bul(60) or b + e 2 G bul(60), and for each 
b G bul(40), either b G bul(60) or b + e 3 G bul(60). Moreover, each bullet for 60 is the 
“image” of a bullet for 54, 51 or 40 in this way. The resemblance to Lemma 3.1 here 


is not a coincidence; Theorem 6.2 makes this similarity precise 


Theorem 4.12. Let M be a monoid. Fix an element x G q (M), an irreducible u G M, 
and a bullet u± + • • • + u r G bul(x). 

(i) If U\ + ■ ■ ■ + u r — (u + x) G M, then u\ + • • • + u r G bul(w + x). 

(ii) If U\ + • • • + u r — (u + x) ^ M, then u + Ui + • • • + u r G bul(u + x). 


Proof. Suppose U\ + • • • + u r — (u + x) G M. For each i < r, we also have 

Ui + ■ ■ ■ + U r — (Ui + U + x) = (ui + ■ ■ ■ + U r — X — Ui) — u ^ M 

since u± + • • • + u r — (u* + x) ^ M . This means U\ + • • • + u r G bul(u + x). 

Next, suppose U\ H-f u r — {u + x) f M. We now verify the necessary conditions. 

(i) Since U\-\ -h u r G bul(x), we have u+U\-\ - \-u r —(u+x) = u±+- ■ -+u r —x G M. 

(ii) For each i < r, u + u\ + • • • + u r — (u^ + u + x) = u\ + • • • + u r — {ui — x) ^ M. 

Thus, u + U\ + • • • + u r G bul(w + t). □ 


Definition 4.13. Fix a monoid M, and an irreducible u G M. For x G q (M), the 
u-cover map -0m : bul(x) —>■ bul(u + x) is given by 


0 m (^1 + ' ' ' + u r ) 


Ui + ■ ■ • + u r if Ui + • • • + u r — (u + x) G M 
u + U\ + • • • + u r otherwise 


for each U\ 4-1- u r G bul(a;). When there is no confusion, we often omit the subscript. 















12 


THOMAS BARRON, CHRISTOPHER O’NEILL, AND ROBERTO PELAYO 


Remark 4.14. The cover map was also defined in fT5j in the context of numerical 
monoids, though its domain was restricted to bullets in which u appears. Both the 
extended domain of 0]^ in Definition 4.13 and the extended domain of oom in Defini¬ 
tion T2 are crucial for the results that follow. In fact, Theorem |4. 15 and Corollary |4 .16 
only held for “sufficiently large” elements of M prior to making these extentions. These 
stronger statements are our primary motivation for extending u>m to q (M), as they en¬ 
sure the correctness of Algorithm 5.6 for all (numerical) monoid elements. 


Theorem 4.15. Fix a monoid M , x G q(M), and a bullet U\ + • — \-u r G bul(a:). For 
each j < r, u\ + • • • + u r — Uj G bill (x — Uj ). 

Proof. Since u± + • • • + u r G bul(x), we have 

(i) {ui + • • • + u r — Uj) — (x — Uj) = Ui + ■ ■ ■ + u r — x G M, and 

(ii) {u\ 4-h u r — Uj) — (ui + x — Ui) = u\ H-1 -u r — (ui + x) ^ M for each i ^ j. 

This means U\ + • • • + u r — Uj G bul(a: — Uj). □ 


Corollary |4.16| follows directly from Theorems |4.12| and |4.15[ and will serve as the 
inductive step of Algorithm |5.6| It also justifies use of the term “cover map”. 


Corollary 4.16. If M is a monoid, x G q (M), and u± + • • • + u r G bul(x), then 

bul(rc) = |j0^(bul(a; - u^). 

i<r 


We conclude the section with Proposition |4.17[ which drastically simplifies the base 
case for Algorithm 5.6 (see Remark 5.1 for more detail). 


Proposition 4.17. Fix a monoid M. For each x G q (M), the following are equivalent: 


(i) oj(x) = 0. 

(ii) bul(x) = {0}. 

(iii) 0 G bul(T). 

(iv) —x G M. 


Proof. We show (i) => (ii) (iii) (i) and (iii) (iv). 

(i) (ii): If co(x) = 0, then any bullet for x has length at most 0, so bul(a;) = {0}. 

(ii) (iii): This is clear. 

(hi) 


i 


This follows from Lemma 4.10 


(iii) (iv): This follows directly from Definition 4.2 


□ 


Example 4.18. Let S = (6,9,20) denote the numerical monoid in Example 4.5 
Proposition 4.17 implies uJs(n) = 0 for every n < —43 = —F{S)\ see Figure [2j 
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5. A DYNAMIC PROGRAMMING ALGORITHM FOR CU-PRIMALITY 


In this section, we apply the results of Section [4] in the setting of numerical monoids 
to obtain Algorithm |5.6| for dynamically computing u- values, with significant runtime 
improvements over existing algorithms. Since any numerical monoid S has finite com¬ 
plement in N, we have q (S') = Z and write Z in place of q(S) in what follows. 


Remark 5.1. One of the main difficulties in developing Algorithm |5.6| was determining 
where to start the inductive process. While Algorithms |3.3| and |3.7| naturally begin at 
the irreducible elements of S, the cu-value and bullet set of an irreducible element x 
are only trivial when x is prime, and this is never the case for non-trivial numerical 
monoids (see Proposition 5.2). One possible (though inelegant) solution is to use 
existing algorithms to compute the bullet sets of every element below some threshold, 
and then use a dynamic algorithm for all further values. 

Thankfully, the extended notion of (u-primality presented in Section [4] allows us to 
avoid this caveat entirely. Proposition |5.5^ t states that for any numerical monoid S, the 
smallest element of q (S') = Z on which ujg takes a nonzero value is —F(S), providing 
a natural base case for Algorithm |5.6[ 


As a non-trivial numerical monoid S contains no prime elements, cug only takes on 
values strictly greater than 1 on elements in S. Proposition |5.2[ however, characterizes 
which elements of q (S) have an cug-value of 1. See [9j, ITS), 1?J for more background on 
pseudo-Frobenius numbers. 


Proposition 5.2. Fix a numerical monoid S = (n i, ..., rik). For eachn G Z, following 
are equivalent. 

(i) w(ra) = 1. 

(ii) bul(ra) = (ei,... ,e fc }. 

(iii) —n is a pseudo-Frobenius number of S. 


Proof. We show (i) 44- (ii) and (ii) <=> (iii). 

(i) =>- (ii): Suppose oj{n) = 1. For each i G (1,we have G bul(n) for 
some bi > 0 since numerical monoids are Archimedian, and by Proposition 4.17[ bi > 0. 
Since u>(n) = 1, each b t — 1. 

(i): This follows from Proposition |4.9 


(ii) =>• (iii): Suppose bul(n) = {e!,...,e fc }. For each i G {l,...,fc}, we have 
Wj-n G S' and n, — n — n, = — n (f S, meaning — n is a pseudo-Frobenius number of S. 

(iii) (ii): If —n is a psuedo-Frobenius number of S, then — n (f S and ?7.j — n G S 
for each % G (1,..., k}, so e* G bul(n). Lemma 4.10 ensures no other bullets exist. □ 

Remark 5.3. Bullets in numerical monoids are usually denoted either as tuples or as 
a list of irreducibles, the former of which is well-suited for storage when implementing 
algorithms. However, for the purposes of implementing Algorithm 5.6, each bullet 
a = (di,..., ak) G bill (n) can be represented using only its value v = a\Ui H-b a^k 
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and length £ = aH-ba*, (a dynamic bullet , see Definition 5.4). Indeed, v is sufficient 

for determining the image of a under cover maps, and £ is sufficient for computing c o(n) 
once all bullets have been found. This is the content of Proposition |5.5| 

Since several distinct bullets can be represented by the same dynamic bullet, the 
cardinality of the computed bullet set with this representation is significantly reduced. 
Much like dynamically computing length sets instead of factorization sets greatly im¬ 
proves efficiency when computating delta sets (Remark |3.10[ ), using this compact bullet 
representation significantly reduces the runtime and memory footprint when dynami¬ 
cally computing w-values. 


Proposition |5.5| is the analog of Lemmas 3T and |3.5| for Algorithm A6 First and 


foremost, parts (b) and (c) prove the correctness of Algorithm 5.6| by ensuring that dy¬ 
namic bullets (Definition |5.4[ ) are sufficient for computing values. Moreover, part (a) 
provides a natural starting place for the inductive procedure in Algorithm 5.6 Lastly, 


part (d) proves that we may further restrict our attention to maximal dynamic bullets, 
and Corollary 6.6 demonstrates the benefit of this reduction (see Remark 6.5). 


Definition 5.4. Fix a numerical monoid S = (n±,... ,nk) and an element n G q(S). 
A dynamic bullet for n is an ordered pair (v,£) G N 2 such that v = Yli=\ biJii and 
£ = are the value and length of some b G bul (n) , respectively. The set of 

dynamic bullets of n is denoted buF(n). The value and length of a dynamic bullet 
(v,£) G buF(n) are the values v and £, respectively, and (v,£) is maximal if £ is 
maximal among dynamic bullets with value v. 


Proposition 5.5. Fix a numerical monoid S = (rii,. .. ,rik), and fix i G {1,..., k} and 
n G q(S). Let if : bul(n) —> bul(n+nf) denote the n*- cover map, : bul(n) —> bnP(n) 
denote the map given by 

b G bul(n) i—> b i n ii Yli =i b i ) e buF(n) 

and ijj* : buP(n) —> bul* (n + nf) denote the map given by 


(v, £) G buF(n) 


(v, £) if v — (rii + n) G S 

(v + rii, £ + 1) if v — (fii + n) f: S. 


(a) If n < —F(S), then bul(n) = {0}. 

(b) For all n G q(5), uj{n) = max{u ! : (v,w) G bul*(n)}. 

(c) For all n G q (S') and i G {!,..., k}, the following diagram commutes: 


bul(n) - > bul(n + nf) 


<f>n+ni 


bul* (n) 


■> bul*(n + nf) 
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(d) A dynamic bullet (v,£) is maximal if and only ifif*(v,£) is maximal. 


Proof. First, Theorem 4.12 ensures that the map if* is well-defined upon observing 
that the value and length of the image of a bullet b under if depend only on the value 
and length of b. Now, parts (a) and (b) follow directly from Propositions 4.17| and 4.9 


respectively, and part (c) follows immediately from Definitions 4.13 and 5.4 Lastly, 


part (d) follows from part (c) and the observation that applying if* to dynamic bullets 
(v,£i), (v,£ 2 ) G bul*(n) preserves any ordering on the second component. □ 

Algorithm 5.6. Given n G S = (ni, ..., nk), computes ca(m) for all m G [0, n] fl S. 
function OmegaPrimalityUpToElement(S', n) 

B rn 4— {(0, 0)} for all m < —F(S ) 
for all m G {—F(S ),..., n} do 
B m ■G- {} 

for all i G {1, 2,..., k} do 

for all (v, £) G B m _ n . do 

{(u, £)} if v — (ni + m) G S 

{(u + Hi, £ + 1)} otherwise 

end for 
end for 

B m {(u,£) £ B m : (v,£) is maximal} 

if m G S then 

co[m] G- maxjf : (v,£) G B rn } 

end if 
end for 
return co 
end function 


Bm t B m U 


Remark 5.7. Table [3] gives a runtime comparison between Algorithm 5.6 and [31 
Proposition 3.3], the algorithm previously implemented in the GAP package numericalsgps. 


Not only is Algorithm |5.6| significantly faster, it also computes an extensive list of co- 
values, rather than just a single value. 


6. Computing w-primality from factorizations 


There are many similarities between Algorithms T3 and A6 beyond their dynamic 
nature. In non-precise terms, bullet sets and factorization sets seem to behave in a very 
similar manner. In this section, we make this connection explicit for numerical monoids 
S by characterizing the bullet set bul(n) of any element n G S in terms of factorization 
sets of certain submonoids of S. As an application, we greatly improve the existing 
bound [n>] on the start of quasilinear behavior (Theorem 6.8) of the cu-function on S'; 
see Theorem 16.121 
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6.1. Runtime of Algorithm 5.6 , We first provide the definition of an Apery set, 
which will be utilized in Theorem 6.2 below. See pj3] for a full treatment of Apery sets. 


Definition 6.1. Let M be a monoid. The Apery set of x G M is defined as 

Ap(M, x) = {m G M : m — x G q (M) \ M}. 

Theorem 6.2. Fix a reduced finitely generated monoid M, and let G = {u\,... ,Uk} 
denote the set of irreducible elements of M. For A C G nonempty and x G M, write 

Z a(x) = {a G Z(x) : ai = 0 whenever Ui ^ A} C Z(x), 

that is, the factorizations in Z M (x) corresponding to factorizations in Z(a)(x). Then 

bul(x) = {be Z A (y + x) : 0 A C G and y G f] Ui£A Ap(M;Ui)} 

for all x G q (M) with ui(x) > 0. 

Proof. First, fix A C G nonempty, y G Pl^eA Ap(M; uf), and b G Z A {y + x). It follows 
that Y^j =1 bj u j — x = y G M, and 

k 

bjUj — (ui + x) = y — Ui M 
3 = 1 

for each Ui G A, so b G bul(x). 

Now, fix b G bul(x). Let A = {ui : 6, ^ 0} C G and y = Y^ k j=i^j u j ~ x ■ Since 


u(x) > 0, A is nonempty by Proposition 4.17 Rearranging the equation for y gives 


y + x = bjUj , meaning b G Z A (y + x), and since y G M and y — Ui^M for 6 * 7 ^ 0, 
we have y G C\ Ui eA Ap(M; uf). □ 

Example 6.3. Let S = (6,9,20) and n = 60 G S'. Our goal is to compute bul(60) 
using Theorem 6.2 For any subset A C { 6 , 9, 20} such that gcd(A) > 1 and 60 G (A), 
each m G f\eA Ap(*S'; a) is relatively prime to gcd(A), so Z A (G0 + m) — 0. The 
only subsets not satisfying both of these conditions are A = {9} and A = {6,9,20}. 
Since Ap(S'; 6 ) D Ap(5'; 9) D Ap(S'; 20) = {0} and the only nonzero m G Ap(S'; 9) with 
Z{ 9}(60 + m ) nonempty is m — 12 , we compute bul(60) as follows: 

bul(60) = Z {6)9] 2 0} (60 + 0) U Z {9} (60 + 12) 

= {(4,4,0), (7, 2, 0), (10, 0, 0), (1, 6 , 0), (0, 0,3)} U {(0, 8 ,0)}. 


Remark 6.4. For each element n in a numerical monoid S, Theorem 622 gives a way 
to compute the set bul(n) (and thus w>(n)) by computing factorizations of elements 
in certain submonoids of S. Since factorization sets in numerical monoids can be 
computed rather quickly in GAP |8j, this yields an algorithm for computing (u-primality. 
This turns out to be faster than [3} Proposition 3.3], the algorithm widely used to 


compute (u-primality in numerical monoids prior to Algorithm 5.6, but not as fast as 


Algorithm 5.6, as it still relies on computing bul(n), rather than the (substantially 
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s 

n 

uj{n) 

[3] Prop. 3.3] 

Algorithm 5.6 

(6,9,20) 

1000 

170 

61319 ms 

6 ms 

(11,13,15) 

1000 

97 

10732 ms 

5 ms 

(11,13,15) 

3000 

279 

874799 ms 

15 ms 

(11,13,15) 

10000 

915 


42 ms 

(15,27, 32,35) 

1000 

69 

234738 ms 

9 ms 

(10,12,15,16,17) 

500 

52 

14338290 ms 

40 ms 

(10,12,15,16,17) 

50000 

5002 


4084 ms 

(100,121,142,163,284) 

25715 

308 


27449 ms 

(1001,1211,1421,1631,2841) 

357362 

405 


3441508 ms 


TABLE 3. Runtime comparison for computing the uj -value of n £ q(S'). 
All computations performed using the GAP package numericalsgps [8] 
and do not make use of eventual quasilinearity. 


smaller) set bill* (n). Additionally, algorithmic use of Theorem 6.2 only produces a 


single cn-value, whereas Algorithm 5.6 produces an exhaustive list of (u-values. 


Remark 6.5. One immediate consequence of Theorem 6.2 is that the number of dis 


tinct values occuring for dynamic bullets in bul*(n) is at most | U t= iAp(S;7u)|. In 


particular, this bound does not depend on n. Since Algorithm 5.6 only stores a single 


dynamic bullet with each value (Proposition 5.51), the inductive step runs in constant 


time in n (for fixed S). This yields Corollary 6.6 


Corollary 6.6. For fixed S, Algorithm F6 runs in linear time in n. 


6.2. Bounding the eventual quasilinearity of (u-primality. For numerical monoids, 
the u;-function admits a predictable behavior for sufficiently large values, a result that 
appeared as [T5| Theorem 3.6] and independently as [TT[ Corollary 20]. Before stating 
this result here as Theorem |6.8[ we give a definition. 

Definition 6.7. A map / : N —> N is quasilinear if f(n) = a\(n)n + ao(ri), where 
ai, a 2 : N —> Q are each periodic and ai(n) is not identically zero. 


Theorem 6.8. [151 Theorem 3.6] Fix a numerical monoid S = (ni,... ,nk). The u- 
primality function is eventually quasilinear. In particular, there is a periodic function 
a(n ) with period n\ and integer N 0 such that for n > N 0 , 

us{n) = —n + a(n ). 

ni 


Remark 6.9. The value of N 0 described in Theorem 6.8 which gives an upper bound 
for the start of the quasilinear behavior of w, is explicitly given in [15]. The actual start 
of the quasilinear behavior, which is called the dissonance , is very often significantly 
lower. Theorem 6.12 will greatly improve the bound for the dissonance. 
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FIGURE 3. A plot showing the function ujs for the numerical monoid 


S = (4,13,19) (left) and (6,9,20) (right), discussed in Example 6.10 


Example 6.10. Figure [3] plots the oj- values for the numerical monoids (4,13,19) and 
(6,9,20). The eventual quasilinearity ensured by Theorem 6.8 manifests graphically 
as a collection of discrete lines with identical slopes, and is readily visible from these 
plots. The proof of Theorem 6.8 appearing in [15j gives a precise (but computationally 


impractical) bound on the start of this quasilinear behavior, and Theorem 6.12 below 
drastically improves this bound; see Table [4] for a detailed comparison. 

We now apply Theorem |6.2 to improve the bound on the dissonance point of the oj- 
function. As with the original bound on the dissonance point given in [15] , this result is 


heavily motivated by computational evidence, in this case provided by Algorithm 5.6 
Lemma 6.11. Fix a numerical monoid S = (ni,n 2 ,... , 14 ). For each n G S, we have 

M(n) < n/n\ < oj(n), 

where M(n ) is the maximum length of a factorization for n. 

Proof. Fix 61 > 0 such that b = Jqei G bul(n) (note that b\ exists since numerical 
monoids are Archimedian). This gives uj{n) > |b| > n/n\. Additionally, we have 
n = Xq=i a j n j — a i n i f° r each factorization a G Z (n), meaning M{ri) < n/n\. □ 

Theorem 6.12. Fix a numerical monoid S = ( n \,..., n ^), and suppose n > N 0 , where 

_ F(S ) + n 2 

° n 2 /ni - 1 ' 

Any maximal bullet b G bul(n) satisfies b± > 0, and uj{n ) = c o{n — n\ ) + 1. 

Proof. Fix a bullet b G bul(n), and let A = {rij : 6, > 0}. First, suppose = 0. 

we have b G Z [n + s ) for some s G p| gA Ap(S'; nf). I 11 particular, 


By Theorem 


6.2 


mm 


(A) f. S, so s — rnin(A) < F(S). Since n > N 0 , we have 


| b | < M(A)(n + s) < 


n + s n + F(S ) 


< 


min(A) min (A) 


+ 1 < 


n + F(S) 


ri2 


+ 1 < n/ni < u>s(n). 
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s 

Dissonance 

Theorem 

6.12 

n 0 pJ)J Theorem 3.6] 

(6,9,20) 

12 

104 

37,800 

(10,12,15) 

190 

325 

66,600 

(10,12,15,16,17) 

10 

175 

34,272,000 

(10,12,13,14,15,16,17,18,19,21) 

10 

115 

450,087,321,600 

(100,121,142,163,284) 

100 

25715 

64,426,520,664,000 

(1001,1211,1421,1631,2841) 

1001 

357362 

« 6 • 10 19 


TABLE 4. Comparison for the start of quasilinear behavior of the u- 
function. The dissonance points were computed using Algorithm |5.6 


available in the next release of the GAP package numericalsgps [Sj. 


by Lemma 6.11 , meaning b is not maximal. This proves the first statement. 

Now, suppose b is maximal, and let a = b — The above argument implies b\ > 0, 


so by Theorem 4.15, we have a € bul(n — nf). This gives 

uj{n) — 1 > u(n — ni) > |a| = |b| — 1 = uj{n) 


1 , 


where the hrst inequality follows from Theorem 4.12 ii) 


□ 


Remark 6.13. The bound given in [15] for the dissonance point of the cu-function 
of a numerical monoid S, while explicit, is far too large to be feasable reached in 


computation. In contrast, the new bound N 0 given in Theorem 6.12 can actually be 
reached computationally; that is, it is often possible to compute every a;-value below 
N 0 (using Algorithm 5.6, for instance). See Table [4] for a detailed comparison of these 
bounds. The a;-value of any element above Nq is determined by the cu-values of the 
elements just below No, and can be obtained by simply evaluating a linear function. 
This makes it possible to compute the cu-value of any element of S. 


7. Future Work 

It is natural to ask whether the dynamic algorithms given above can be extended 
to settings outside the realm of numerical monoids. In fact, Algorithms 3.3 and 3.7 
extend immediately to any finitely generated monoid M, once a suitably “bounded” 
subset of M is chosen (for instance, one co uld use the set of divisors of some element 
m G M). In order to generalize Algorithm 3.9, the eventual behavior of A : M —>■ 2 1 
must be examined. 

Problem 7.1. Describe the eventual behavior of A : M —> for M finitely generated. 


It is worth noting that the generalization of Algorithm 3.7 to finitely generated 


monoids does allow the use of computation packages for investigating Problem 7.1 


Generalizing Algorithm 5.6 presents a slightly more subtle issue. Section [4] extends 
the ^-function on a finitely generated monoid M to its quotient group q(M), along 
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FIGURE 4. The monoid M CN 2 given in Example F2 The elements of 
M are marked by filled dots, the lightly shaded region marks the divisors 
of m = (3,1), and the darkly region marks —M. 


with both the 
an element m 


see Example 7.2 


“base case” and the “inductive step” for Algorithm A6. However, given 
G M, the set of divisors of m with nonzero cn-value need not be finite; 
Thus, more care must be taken in order to dynamically compute 


(n-values in M, and it is as yet unclear how to make this distinction. 

Example 7.2. Consider the monoid M cN 2 consisting of all points not lying above 
the ray generated by (2,1). Definition 4.2 extends the domain of ojm to all of Z 2 , and 


Proposition 4.17 states that = 0 precisely when m G —M. However, the set of 

divisors of any nonzero element m G M contains infinitely many elements of Z 2 which 
necessarily have both nonzero oj -value and nontrivial bullet sets; see Figure [4} 


Problem 7.3. Develop a dynamic programming algorithm to compute u-values for 
finitely generated monoids. 

Given an element n G S, its catenary degree c[n ) can be computed from a graph 
with vertex set Z (n) j6]. Like delta sets, the catenary degree function c : S —» N is 
periodic for large input values [5], but no bound on the start of this periodic behavior 
is known, and the period is only known to divide the product of the generators of S. 
The iterative nature of the periodicity proof given in [5] hints at the existence of a 
dynamic algorithm for computing catenary degrees in numerical monoids, allowing for 
the use of computation when investigating the eventual periodic behavior. 

Problem 7.4. Find a dynamic programming algorithm to compute the catenary degrees 
of numerical monoid elements. 
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